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ABSTRACT 

The Kepler Mission recently identified 997 systems hosting candidate extra- 
solar planets, many of which are super-Earths. Realizing these planetary systems 
are candidates to host extrasolar asteroid belts, we use mid-infrared data from 
the Wide-field Infrared Survey Explorer (WISE) to search for emission from dust 
in these systems. We find excesses around eight stars, indicating the presence of 
warm to hot dust (~100-500 K), corresponding to orbital distances of 0.1-10 AU 
for these solar-type stars. The strongest detection, KOI 1099, demands ~500 K 
dust interior to the orbit of its exoplanet candidate. One star, KOI 904, may 
host very hot dust (~1200 K, corresponding to 0.02 AU). Although the fraction 
of these exoplanet-bearing stars with detectable warm excesses (~3%) is similar 
to that found by Spitzer surveys of solar-type field stars, the excesses detectable 
in the WISE data have much higher fractional luminosities (Ld us t/L*) than most 
known debris disks, implying that the fraction with debris disks of comparable 
luminosity may actually be significantly higher. It is difficult to explain the pres- 
ence of dust so close to the host stars, generally corresponding to dust rings at 
radii <0.3 AU; both the collisional and Poynting-Robertson drag timescales to 
remove dust from the system are hundreds of years or less at these distances. 
Assuming a steady-state for these systems implies large mass consumption rates 
with these short removal timescales, meaning that the dust production mecha- 
nism in these systems must almost certainly be episodic in nature. 



Introduction 



The evolution of solid circumstellar material is a competition between the accretional 
processes that form planets (when speeds are low) and the violent collisional processes which 
in post-formation systems inexorably grind down the material into particles small enough 
that radiation forces can remove them from the system. The disappearance of infrared 
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(IR) excesses around forming stars (on 3-10 Myr time scales) occurs because the accretional 
processes win and sequester the huge mass of dust into a (relatively) tiny number of objects 
with very low surface area/mass ratios. In contrast, the observed debris disks (known around 
hundreds of solar-type stars) are thought to be from collisional dust production as velocities 
rise during the final stages of planet assembly (for relatively young stars), or triggered around 
older stars as dust is liberated in higher-speed collisions (jWyattl 120081 ) . It is these more 
mature systems where our understanding of the physical picture is arguably incomplete. 
Althoug h only ~2% of nea rby main-sequence stars have massive detectable warm disks 
(>200K; lLawler et al.l 120091 ). several of these disks require dramatic hypotheses to produce 
the large dust mass estimated. That is, estimates of the collisional grinding of small-body 
belts indicate that the systems in qu estion sh ould not sti l l poss ess dust at anywhere near 
the observed level ( jWyatt et al.ll2007l ). In fact, iRhee et al.l ( 120081 ) went so far as to liken the 
dust mass around BD +20 307 as a 'miracle' given collisional-model predictions. 

Debris disks are detected and studied via the excess IR flux they produce in the system, 
caused by the dust re-radiating the absorbed stellar light. This dust emission can then 
dominate the stellar photosphere in the IR, where the dust's blackbody peak is set by the 
equilibrium temperature at the orbital distance (for a narrow ring) or at an extended disk's 
inner edge, because the warmest dust usually dominates the systems dust emission. For 
'realistic' (although uncertain) dust grain size distributions, cross- sectional area and thus 
thermal emission is dominated by small grains, but the system mass is dominated by the 
unconstrained larger bodies. The true observable is thus the mass in grains in the decade 
larger than the peak wavelength. 

Known disk systems consist of two types. The most common are massive and cold, 
analogous to the Solar System's Kuiper belt. A few hotter (T>200K) systems are known, 



which are massive cousins of our asteroid belt. An examp 



e is HP 69830, hosting 3 RV- 
2006h and a wa rm dust ring at 



discovered Neptune-mass planets inside 0.6 AU ( lLovis et al.l 

~1 AU, with an estimated 4xl0~ 7 M ffi in small grains alone (jLisse et al. 2007 ). Production of 



this dust is problematic ( IBeichman et al.ll2005l ). with hypotheses ranging from comet swarms 
to super-comets to massive planetary collisions; it is difficult to judge the plausibility of these 
mechanisms given the rarity of these outcomes. 

Dust grains in orbit around a star are quickly destroyed, and must be replenished by 
a fragmentatio n cascade, wh ere asteroids within a belt collide and break into progressively 
smaller pieces ( jWyattl 120081 ). The lifetime of a debris disk will be the same as the lifetime 
of the largest bodies in the source population that are available to be pulverized into dust 
grains. Very small grains (<1 /im) are quickly blown away by stellar radiation pressure. 
Larger grains suffer Poynting-Robertson (PR) drag and spiral into the star in much less 
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than the system's age; for dust very close to the star (~0.1 AU), orbits collapse in only 
hundreds of years. 

In this paper we search for excess infrared emission, indicating the presence of a dusty 
debris disk, in systems found by the Kepler mission to possess transiting exoplanets. In 
Section [2] we define our sample and discuss our methodology for identifying excess candidates, 
and in Section [3] we provide additional details on the eight stars that have excesses. Section H] 
contains a discussion of the problems associated with having this much dust so close to 
the host star, and finally Section contains a summary of recent work and this study's 
contribution relating debris disks to the presence of exoplanets. 



2. WISE Data 

WISE released preliminary survey data in April 2011, covering about half the sky in 
photometric bands centered on 3.4, 4.6, 12, and 22 /zm (Wl, W2, W3, and W4, respectively; 



Wright et al.ll2010f ). Due to the pointing strategy of the telescope, the coverage on the sky is 
non-uniform. Many individual images are stacked to produce the preliminary survey data, 
and some areas of the sky were visited more often than others, providing deeper coverage in 
those areas. 



2.1. Defining Our Sample 



We use WISE data to search for dust emission around the 928 solar-type stars (F, G, and 
KO-7 spectral types; T e fj=4000-7200 K) that have stellar data (TW, metallicity, and log(g)) 
reported in the Kepler database as of November 2011 and are candidates to host one or 
more p lanets contained in the first Kepler data release (IBorucki et al.ll2011t iLissauer et al. 
201 lbl ). Although most of these transiting planet candidates remain to be independently 
confirmed, we will refer t o them as planets for brev ity. The false-positive rate for single 



planet s ystems is <10% (IMorton and Johnson! 12011 ) , and is much lower for multiplanet 



systems (IRagozzine and Holman 



201 



0: 



Lissauer et al. 



2011aJ). In any case, the existence or 



lack of planets is irrelevant for our analysis of IR excess. 

There are several M stars reported as having planets in the Kepler database, but we 
did not consider these for analy sis because of known systematic problems w ith modeling the 
atmospheres of such cool stars (jBeichman et al.ll2006al ; ISinclair et al.ll2010l ). We also did not 
consider the few hotter (A-type) stars with known Kepler planets in order to keep to stars 
similar to our Sun. 
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The WISE preliminary data release covers about half the Kepler field, detecting 439 
out of the 928 solar-type stars; we consider a target to be 'detected' in a given band if the 
WISE Source Catalog reports a signal-to-noise ratio (SNR) greater than 3. (WISE does not 
report fluxes, but rather magnitudes and their corresponding magnitude uncertainty.) All 
WISE detections are within 1 arcsecond of the coordinates given by the Kepler team. The 
band wavelengths cover where hot dust emission would peak, which in these systems is also 
near the equilibrium temperatures of the known Kepler planets, most of which are located 



0.01-0.4 AU from their host stars (ILissauer et al.ll2011bl ). 



134 of these stars were flagged by the WISE team as having contamination from nearby, 
bright stars in the form of a halo or diffraction spike, or flagged as extended, suggesting 
probable contamination by a background object. These stars were not considered further. 

The remaining 325 stars are all detected in both WISE's Wl (3.4 /j,m) and W2 (4.6 /zm) 
bands, ~2/3 of these are also are detected in W3 (12 jttm), and only six are also detected in 
W4 (22 fjm). 



2.2. Finding Excesses in the WISE Data 



To determine if the WISE data match the expected stellar photosph eric spectrum or rep 
resen t an excess of emission, we use Kurucz stellar atmosphere models (ICastelli and Kurucz 
20041 ) chosen from a grid of stellar temperatures, surface gravities, and metallicities to closely 
match those given for each star in the Kepler online database (IBrown et al.ll201ll ). The 
Kurucz models are s caled using photometry from the Sloan Digit al Sky Survey (SPSS ; 
Abazaiian et alil2009h and the Two Micron All-Sky Survey (2MASS; Iskrutskie etaDbood ). 
These Kepler targets are more distant than most known debris disk hosts (hundreds of par- 
sees), so extinction must be taken into account. The Kepler database provides E{B — V) 
meas urements for each star, which we scale for each photometric band ( IRieke and Lebofsky 
19851 ) and use to correct the photometry The Kurucz model is scaled to each star's flux 
using the corrected photometry, and then convolved with the efficiency function over each 
WISE passband. The band zero points then provide expected magnitudes. 

There are several sources of error in these measurements. To estimate magnitude un- 
certainties 5, the given error bars for each datapoint (5w) and the WISE abso lute photo- 
metric errors (2.4, 2.8, 4.5, and 5.7% in Wl-4, respectively; IWright et al.ll2010l ) are added 
in quadrature. To assess the significance of each possible excess, the WISE magnitude (m) 
is subtracted from the predicted stellar magnitudes (m pre d), and this is divided by the un- 
certainty (5) in that band to give a dimensionless excess. Our sign choice results in positive 
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Fig. 1. — The excess significances for each band (m pre ^-m) / 5 . Overplotted are Gaussians 
with width 5 for reference. Small (~1 5) systematics are still clearly present. Filled bars show 
what we accept as significant excesses. One system (KOI 379) shows excess in Wl and W2 
above the 5 S level (7 S in Wl and 5 6 in W2) but is not considered to be significant because 
of an obvious large flux mismatch between the SDSS and 2MASS photometry, perhaps 
indicating stellar variability. 



excesses indicating that the WISE flux is higher than the predicted flux. The uncertainties 
in the stellar measure ments reported by the Kepler database are ±200 K in temperature 



and 0.4 dex in log(g) ( IBrown et al.ll201lf ). We estimated the contribution of these uncer- 
tainties to the significance of the WISE excesses by using stellar models from (T eff , log(g)) 
gridpoints that were varied by the error bars in each of these parameters, and we found that 
this resulted in <1 5 changes in the excess. In order to account for this without having a 
formal number to add to our errors, we adopt the fairly stringent requirement that in order 
for an excess to be considered significant, it must be more than 5 5 above the photospheric 
magnitude. 

Because of the large photometric aperture on the WISE measurements (8.25 arcseconds 
for W1-W3 and 16.5 arcseconds for W4), we visually inspected both the WISE and SDSS 
images for each of the stars we find to have an excess and found that there is not contam- 
ination from any nearby, bright stars. In addition, we checked to make sure there were not 
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Fig. 2. — Histograms of each of the four stellar parameters reported by the Kepler team. 
Filled bars show these values for the eight systems with significant excess in at least one 
WISE band. There appear to be no correlations between the presence of a significant excess 
and any of these parameters. 

systematic offsets in the expected position of the stars with excess, indicating contamination 
by a nearby background source. 

The distribution of excesses for each band are shown in Figure [TJ They follow approx- 
imately Gaussian distributions, with very few datapoints falling above 5 5 and none falling 
below -5 5. Figure [1] makes it clear that the systematic error here actually underestimates the 
significance of the excesses; if anything we are being more conservative with our definition 
of what is a significant excess and what is not. 

Figure [2] shows stellar properties for our sample as reported by the Kepler team: ex- 
tinction, stellar temperature, surface gravity, and metallicity. No trends are seen between 
any of these data and the presence of a significant excess in any band. 

3. Candidate Debris Disk Detections 

The vast majority of stars we examined have WISE fluxes consistent with the photo- 
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Table 1. FGK Stars with No Significant (>5 S) Excess 



KOI 



1 


127 


257 


379 


488 


574 


678 


772 


864 


942 


1094 


1169 


1288 


1403 


1489 


5 
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258 


384 


490 


579 
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775 
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1302 


1404 


1501 


7 
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274 


385 


496 


583 
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1202 
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10 


151 
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392 


503 
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on 

20 


loo 


281 


A no 

408 


506 
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588 
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800 
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no 1 
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1 1 no 

1108 
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1 A fW 

1407 


1515 


22 


176 


283 


409 


508 


589 


701 


801 


884 


984 


1109 


1216 


1311 


1408 


1519 


69 


180 


284 


416 


510 


600 


704 


806 


889 


986 


1110 


1218 


1312 


1409 


1521 


72 


188 


292 


419 


512 


605 


707 


809 


890 


992 


1115 


1219 


1314 


1413 


1525 


75 


189 


296 
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513 
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709 
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891 


998 


1116 


1220 


1337 


1419 


1526 


82 


194 


297 


431 


518 
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712 


822 


898 


999 


1117 


1236 


1338 


1423 


1529 


84 


201 


299 


442 


523 


628 


717 


824 


900 


1001 


1118 


1240 


1339 


1427 


1530 


98 


202 


301 


443 


524 


632 


718 


830 


901 


1010 


1128 


1244 
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1430 


1533 


100 


209 


304 


444 
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644 
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834 


902 


1022 


1129 
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1432 


1534 
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313 
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837 
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1434 
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1367 
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1273 
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1445 


1588 
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227 


345 


472 


537 


661 


755 


853 


928 


1054 


1150 


1276 


1387 


1465 


1605 


118 


232 


349 


474 


554 


662 


759 


855 


929 


1061 


1160 


1278 


1395 


1477 


1608 


122 


241 


351 


475 


557 


663 


765 


858 


934 


1086 


1168 


1281 


1396 


1480 


1609 


123 


242 


374 


479 


561 


674 





















Table 2. Debris Disk Candidates 



Kepler Data 10 BB Dust Ring Model" 1 /nn Mod. BB Dust Ring Model 6 







Stellar 




Planetary 






Dust Properties 






Dust Properties 




KOI 


sp. 


extinction 


a 


radius 


mass 1 * 


T 


a 


M 


Ld/L* 


T 


a 


M 


L d /L, 




type c 


E(B-V) 


[AU] 






[K] 


[AU] 


[M 9 ] 


[xl0~ 3 ] 


[K] 


[AU] 




[xl0~ 3 


904 e 


K5 


0.074 


0.029 


2.1 


5 


1200-f 


0.02 


lxl0~ 8 


13 


noo-f 


0.02 


2xl0~ 9 


15 








0.159 


3.0 


10 


















469 


GO 


0.116 


0.095 


5.5 


35 


400* 


0.36 


8xl0 -7 


3 


3O0* 


0.64 


2xl0 -7 


3 


559 


K0 


0.099 


0.052 


1.4 


2 


471 


0.23 


4xl0~ 7 


4 


382 


0.34 


lxHT 7 


1 


871 


G5 


0.092 


0.105 


10.9 


144 


329 


0.27 


6xl0~ 7 


4 


283 


0.37 


lxl0~ 7 


4 


943 


K0 


0.102 


0.045 


2.2 


5 


357 


0.28 


lxl0~ 6 


6 


303 


0.39 


2xl0~ 7 


6 


1099 


G5 


0.114 


0.573 


3.7 


15 


487 


0.14 


6xl0~ 7 


14 


394 


0.22 


lxKT 7 


14 


1020 


G3 


0.11 


0.294 


21.9 


614 


lOO-f 


9.89 


5xl0 -4 


3 


90 } 


12.21 


6xl0" 5 


2 


1564 


G4 


0.109 


0.275 


3.1 


11 


135-f 


1.96 


lxl0~ 4 


16 


115* 


2.70 


2xl0 -5 


14 



a Assuming blackbody dust grains of radius 10 /im 

b Assuming modified blackbody dust grains of radius 1 (im 

c Spectral type is approximate and based on the stellar temperature provided by the Kepler team 
d Planct mass is calculated according to M/M® = (R/R m ) 2 - 08 l lBorucki et alj|201 it) 
c KOI 904 has two known planets, both of which are listed here 
f Td us t is assumed rather than fit. See text for details on each system. 
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wavelength |>m] 





wavelength [yim] 



Fig. 3. — Nine stars with WISE data that is consistent with the stellar photosphere. Included 
here are all of the confirmed systems that currently (as of November 2 011) have WISE dat a 
(none of these syste ms have significant excesses): KOI 7 = K epler 4 feorucki et al. 2010 ). 
KOI 10 = Ke pler 8 faorucki et al.l fioioh. KOI 20 = Kepler 12 jFortnev et al.lboill) . KOI 72 



Kepler 10 flBatalha et al.ll201ll) . KOI 84 = Kepler 19 



Kepler 7 faorucki et al.lboiol ). and KOI 377 = Kepler 9 Jrlolman et al 



Ballard et al 



2011) . KOI 97 = 
2010l ). KOI 1427 



and KOI 258 are included as additional examples, being one of the coolest and warmest 
stars in the sample, respectively. Blue curve represents the stellar spectrum model. Cyan 
diamonds and squares show the SDSS and 2MASS data respectively used to scale the stellar 
model, and the red triangle shows the Kepler visual magnitude estimate converted to flux. 
The filled symbols show the extinction-corrected data, and the empty symbols show the 
original measured values. Green filled circles show the WISE data, and the error bars are 
often smaller than the datapoints. The excess (in multiples of S) of the observed-predicted 
magnitudes are given for the stellar model above each WISE datapoint. At the top, triangles 
show the planetary thermal emission wavelength peak for reference. 
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sphere (these are listed in Tableland a selection of these stars are shown in Figure |3]). In 
our full sample, we find eight stars with an excess of 5 5 or higher in at least one band. 
Unsurprisingly, these excesses generally appear in the longest two wavelength bands. These 
stars are listed in Table [2J the WISE data for each of these systems are listed in Table El 
and we discuss each of them in Sections I3.1H3.3I below. 

Many simple debris disk models have be e n made using blac kbody dust grains (i.e. 
Beichman et al.l l2006bl ; iHillenbrand et al.ll2008t lLawler et al.l 120091 ) . We adopt 10 /im ra- 
dius blackbody dust grains for an initial model, both for ease of comparison with previous 
literature and for simplicity given the lack of spectral data. These models are not attempting 
to find the highly unconstrained 'best' possible dust model due to the myriad of complica- 
tions inherent in dust emission (composition, radial distribution, emission features, and grain 
size), but rather to prove that the presence of dust provides a physically plausible model 
that approximately reproduces the WISE measurements. A single temperature blackbody 
is scaled to fit the WISE photometry, giving Ldns t/L*, which i s used to calculate the area of 
absorbing dust o"d U st using adust = 47ra 2 Ld U st/L* (jWyattl l2008f ) . When a density and radius 
are assumed for the grains, this can be converted to the mass in dust. We used 3.3 g cm -3 
for silicate dust grains. The temperature changes the wavelength of the peak of the curve. 
The orbital distance of the dust is then calculated using the same assumptions as were used 
to calculate Kepler planet temperatures (albedo of 0.3, uniform surface temperature, and no 
atmospheric effects; iBorucki et al.ll201ll ). 



It is probable that smaller grains are present, and if they dominate in the size distribution 
of dust grains, the dust could be significantly farther out and produce the same fit. We show 
this uncertainty in Table [21 where in addition to the best-fitting 10 /im dust grain model we 
show the best fit for a disk made of 1 /im modified blackbody dust grains, where the emissivity 
of the grains drops proportional to (A/A ) _/3 , with A equal to the grain size and (3 equal to 
1 (IDent et al.l 120001 ) . This means that the spectral energy distrubution resulting from these 
grains is the same as a blackbody for wavelengths shorter than 1 /im, and drops off more 
steeply than a blackbody at larger wavelengths. Smaller grains may be more realistic, as 
resolved disk studies have found larger disks than predicted by fitting a blackbody spectrum 



to the mid-IR excess (i.e. iKrist et al.ll2010l ). However, in the interest of picking a uniform 



model that can be compared with previous Spitzer studies, we use the 10 /im blackbody dust 
model for most of our calculations. 

For our 10 /im dust grain model, the resulting blackbody curve is added to the stellar 
atmosphere model, and this is used to fit the WISE data. For most systems, the temperature 
and surface area of the emitting dust are both free parameters in the fit. Because in some 
systems there is only one datapoint that rises off the photosphere, the temperature cannot 
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be determined uniquely even if an upper limit can be set. These systems are noted in Table [2] 
and discussed below. For all eight systems, the fit for temperature is fairly unconstrained; 
the longest wavelength datapoint is almost always still rising, allowing a large number of 
cooler blackbody curves from more massive dust disks to pass through the datapoints. Thus 
the temperatures we calculate should be thought of as upper limits, as a larger amount of 
somewhat cooler dust could also produce the observed emission. 

Using this simple dust model, the dust masses needed to produce these excesses range 
from — 10~ 8 -10~ 4 Earth masses, orders of magnitudes larger than the dust mass of our solar 
system (~4x 1(T 10 M a 
systems (typically 10~ 8 — 10 





Hahn et al. 


200: 


2), and comparable in mass to known debris disk 


-10- 2 M ffi ; 


Beichman et al. 


2005; 


Brvden et al. 


2006; 


Wvatt et al. 


Hillenbrand et al. 


2008). but with higher fractional luminosities as 



shown in Figure HI 

We find that our eight excess systems can fit into three categories based on dust tem- 
perature (and thus dust location). For much of the literature, >200 K dust is considered 
to be hot dust. Due to the large variation in debris disk temperatures and terminology, 
for this paper we define our debris disk categories as the following: hot (—1000 K), warm 
(~300-500 K), and cool (—100 K). We now discuss each excess system in turn. 



3.1. Warm Dust Systems (-300-500 K) 

Because of its wavelength range, WISE is most effective at detecting dust temperatures 



between 300-500 K, similar to the handful of warm debris disks discussed in IWyatt et al. 



(120071 ) . Figures \5\ and show the five systems that possess excesses consistent with emis- 
sion by warm dust. Not surprisingly, due to the low photospheric flux predicted at long 
wavelengths, none of these systems have W4 detections with SNR>3. The low SNR data- 
points and upper limits that are plotted (but not included in our fit) give tantalizing hints 
of what the upcoming full WISE data release could reveal about the dust locations in these 
systems. The W4 datapoints may reveal that these systems are actually similar to the cool 
dust systems in dust mass and temperature (see Section |3~2]) . Only one of these five systems 
(KOI 1099) requires dust warmer than >300 K (see discussion below); in the other four 
systems we may be seeing the Wein tail of a cooler, more massive dust ring. 

With the currently available WISE data, we find that these disks are plausibly hotter 
than the majority of known debris disk systems, and our best fitting models place dust on 
orbits of a few tenths of an AU. This is perhaps not surprising, as it is similar to the orbital 
radii of the known planets in these systems, so the dust may be co-located with several Earth 
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Dust Temperatures [K] 




Dust Temperatures [K] 



Fig. 4. — The derived masses, fractional luminosities, and temperatures of dust rings in our 
sample (using 10 /xm radius grains ) comp ared to o ther known systems. Red triangles show 



syste ms examined by IWyatt et al.l (120071 1 (see also iBeichman et al.l 12 005. 



2006bj LisseetaL 



20121 ). green squares are from a Spitzer spectroscopic survey ( 



pentagons are from a Spitzer photometric survey at 70 fim ([Hillenbrand et al.l 12008). Fo r 



>awler et al 



20091) , and green 



comparison, the purple diamond shows zodiacal dust in our solar system ( IHahn et al 



20021 ). 



All systems shown are mature, solar-type stars. Blue circles show the dust models presented 
in this paper. Because M grain is proportional to grain size r grain , using smaller grains would 
result in lower total dust masses (see Table [2]). While the masses shown here are similar to 
known systems, because of the warmer temperatures we find for these Kepler systems, the 
resulting fractional luminosities Ld us t/L* are much higher. 
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KOI 1099 

, , T 




wavelength [/im] 

Fig. 5. — KOI 1099, our best candidate for a debris disk due to the strong excess in W2 and 
W3, and a weak excess in Wl. See caption for Figure|3j The excess values (in multiples of the 
total magnitude uncertainty 8) for the stellar atmosphere model predictions (not including 
the dust model) are shown above each WISE datapoint. Additionally, the purple dotted 
curve shows the blackbody spectrum of a thin ring of 10 /im-radius dust grains, and the 
solid red curve shows the star+dust spectrum. The best-fitting dust mass, temperature, and 
orbital distance are also shown. The unfilled WISE datapoint shows SNR<3 data that was 
not used in the fit, but gives a preview of where the upcoming WISE full data release may 
provide significant data. 



masses of planetary material. Most systems in the literature known to host both a debris 
disk and one or more exoplanets have a large separation (tens of AU) between the orbit of 
the planet (s) and the debris disk (IBryden et al.ll2009i ; iDodson-Robinson et al.ll201ll ). One of 
the five systems in our sample with warm excesses has the modeled dust ring interior to the 
known planet (Figure [5]), with four exterior (Figure |6]). All planets are within a few tenths 
of an AU of the modeled debris disk, which may allow dynamical interactions between the 
planet and the disk. 

The G5 star KOI 1099, shown in Figure EJ is our most promising candidate to host a 
debris disk, appearing to have excesses in all WISE bands. The weak excess in Wl combined 
with highly significant excesses in both the W2 and W3 bands demands dust emission of 
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Fig. 6. — Stars with warm excesses, at orbital distances similar to that of the known planets. 
See captions for Figures [3] and [51 Open green triangles show upper limits on the W4 data, 
while open circles show W4 data with SNR<3. 

~500 K using our nominal model. An excellent match to all photometry between 0.5- 
12 jum is provided by a confined dust ring at 0.14 AU, although more complex distributions 
are certainly possible. The 1 /jm modified-blackbody alternate model results in a cooler 
temperature (~400 K) and larger orbital radius for the dust ring (0.22 AU). Either model 
results in a ring is inside the 0.57 AU orbit of the known Neptune-mass planet (candidate), 
and may point to inspiral of dust past the planet from a more distant dust source. The W4 
flux has SNR<3 and so was not used in our fit, but if future WISE data releases confirm the 
~1.5 mJy flux, this points to an even greater dust mass, much of which would be at larger 
distances from the star. Using the low SNR W4 datapoint along with the other WISE bands 
would not allow a single temperature fit, simply indicating that this system contains dust 
that is producing significant emission at more than one orbital distance. 
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The other four warm systems, shown in Figure [61 do not have formally significant Wl 
and W2 excesses. It is possible that these weak Wl and W2 excesses should be ignored, in 
which case the W3 excess is the start of a blackbody of a much more massive but colder disk. 
We evaluate this explicitly for the KOI 943 system, where there is a low-SNR W4 datapoint 
rather than just an upper limit. 

KOI 469 is a GO type star with a super-Neptune at 0.095 AU. Because only the W3 
datapoint is off the photosphere (with only an upper limit in W4), for our best-fitting dust 
ring we fixed the temperature by eye and fit the data by scaling the mass of the ring. This 
gives a dust ring at 0.36 AU, outside the orbit of the known planet. 

KOI 559 is a K0 type star with a ~2 M® planet at 0.05 AU. There is a very weak excess 
in Wl and W2, with a strong excess at W3 and only a W4 upper limit. The best-fitting 
dust ring sits at 0.23 AU, well outside this hot super-Earth's orbit. 

The G5 star KOI 871 has a Saturn-mass planet at 0.1 AU. The WISE data is similar to 
KOI 559, with a very weak excess in Wl and W2, a strong excess in W3, and a W4 upper 
limit. The best fitting dust ring sits at 0.27 AU, again outside the planet's orbit. 

The K0 star KOI 943 follows the usual pattern of weak excesses in Wl and W2, with a 
strong excess in W3. The best-fitting dust ring sits at 0.28 AU, outside the orbit of KOI 943's 
super-Earth at 0.05 AU. The low-SNR W4 datapoint hints at a more massive, more distant 
dust dust ring that will perhaps be confirmed by the upcoming full WISE data release. 
Using the nominal 1.5 mJy W4 flux in the fit pulls down the temperature to about 150 K, 
corresponding to a ring at 1.5 AU with a mass of 10~ 4 M®, similar to the cool dust systems 
described in Section 13.21 

The dust masses required for all of these systems are fairly low: ~10~ 7 -10~ 6 M®, ap- 
proximately equivalent mass to grinding a 10 km asteroid into 10 jum dust grains. 



3.2. Cool Dust Systems (~100 K) 



These systems (Figure [7]) more closely resemble classical debris disk systems that have 
been studied extensively at longer wavelengths by Spitzer, with du s t located at several AU, 
well outside the orbits of the known planets ( iBryden et al.l l2009t iDodson-Robinson et al. 
201ll ). We find two stars in our sample with excesses that are consistent with blackbody peaks 
at the W4 effective wavelength of 22 /mi or longer, and strongly significant W4 detections. 

KOI 1020 is a G3 star, and is the only star which has a strong W4 excess and yet 
no W1-W3 excess. This would require the dust to be quite cool (< 100 K). Due to the 
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Fig. 7. — Stars with cool excesses (~100 K). See captions for Figures |3] and [5j 



unconstrained temperature fit, we fixed Td us t at 100 K, but it could be colder if a more 
massive disk is present. Even the 100 K temperature places the dust at about 10 AU, far 
outside the orbit of its hot Jupiter at 0.3 AU. 

KOI 1564 is a G4 star, and has weak excesses in each of the 3 shortest wavelength 
bands, with a large excess in W4. Due to the large error bars on the W3 and W4 data, there 
are quite a range of temperatures that fit, but our best dust model has an orbital distance 
of 2.0 AU. KOI 1564's known Neptune-mass planet orbits well inside this, at 0.28 AU. The 
future improvement in the error bars provided by the next WISE data release will certainly 
better constrain this belt. 

Both of these stars require high dust masses (~10~ 4 M e ) to produce the large excesses 
observed in W4, similar to the masses calculated for Spi tzer surveys of solar-type stars at 24 
and 70 /im (IBryden et al.ll2006l ; iHillenbrand et al .1120081 ). although KOI 1564's disk may be 
among the hottest known with this mass or larger. 



3.3. Hot Dust (~1000 K) in a Multiplanet System 



The K5 star KOI 904 (Figure [S]) is the only multiplanet system to show an excess, as 
well as the only system with signficant excess in the shortest wavelength band (3.4 /im). This 
alone requires extremely hot dust. Both of the planets are super-Earths, one at 0.03 AU 
and one at 0.16 AU, indicating that there is at least ~15 M ffi of material near the star. 
Having two planets i s significant, because it greatly reduces the possibili ty of a false positive 
exoplanet detection ( jRagozzine and Holmanll2010t iLissauer et al.ll2011ai ). Currently, there is 
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Fig. 8. — KOI 904: the only star with a strong excess in both Wl and W2, and the only 
multiplanet system in this sample to show an excess. See captions for Figure [31 [5] and |6j 

not a detection in the W3 or W4 bands, only a low-SNR W3 point and a W4 upper limit, so 
the dust model is not well constrained. For lack of any other guidance on the temperature, 
we fixed it so the blackbody curve peaks near the W2 effective wavelength and found an 
acceptable fit. This results in very hot dust, close to the sublimation temperature of silicates. 
This exact temperature is uncertain, but must be quite hot (>1000 K) in order to fit both 
the Wl and W2 datapoints. But with currently only two datapoints, more detailed dust 
modelling is unwarranted. 

This remains our most insecure excess; there is a small mismatch between the SDSS and 
2MASS photometry which could be the result of stellar variability (although the star is not 
flagged as variable in the WISE database). This is among the cooler stars in the sample, and 
so there is the possibility of poor stellar atmosphere modeling, as for the M stars. However, 
other K5 and later stars appear in our sample without showing signs of excess in the WISE 
data (see Figure [2]). The very hot dust is also a bit worrying, though it is not very different 
from the equilibrium temperature of the innermost known planet in this system (960 K), 
and dust temper atures in the range 500-1200 K have also been detected around white dwarfs 



f lChu et aUEOllt ). 
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If real, this very hot dust cannot be produced by planetesimals in orbit at this distance; 
as discussed in more detail in Section HI both the collisional timescale and the Poynting- 
Robertson drag time guarantee that particles cannot stay in such close orbits to the star for 
more than hundreds of years. 

Only a very small dust mass (~10 -8 M©) is required in this system because of the 
extremely hot temperatures. Dust this hot has only been observed with optical interfer - 
ometers, and may actually be fairly common around solar-type stars (IDefrere et al.ll20!ll ). 
but is difficult to observe. Excesses this hot may have been missed by Spitzer spectroscopic 
surveys, which assumed no excess was present at short er wavelengths if no slope other than 



Rayle igh- Jeans was found in the short wavelength data (jBeichman et al.ll2 006b 



20091 ) . This excess is perhaps similar to that found for HD 23514 (jRhee et al 



Lawler et al. 



2008), which 



was also detected at short wavelengths only. The full WISE data release should provide a 
high SNR W3 datapoint; if it remains above the photosphere, the case for a tenuous, but 
very hot ring of small grains near the planets will be strengthened. 

We note that as this paper went to press, the February 2012 Kepler data release 
(jBatalha et al.ll2012l ) provided 3 additional planets in this system, bringing the total number 
of planet candidates around KOI 904 to five, all within 0.2 AU. 



4. Dust Production 

For the systems with abundant dust interior to 1 AU, it is difficult to produce viable 
steady-state models with such a large mass of small grains remaining so close to their stars. 
The PR drag timescale for dust to spiral into the star is extremely short at these close 
distances: 

'"grain \ / 0> \' 



tpR = 200 yr — — 
K J \ 10 urn) VI AU/ 



for 3 g cm -3 grains (IBurns et al.lll979f ). giving only ~10 years for 10 /j,m grains at 0.25 AU 
(a typical orbital distance in Table [2]). Inspiraling grains would not ultimately reach the 
star in any case, as they should vaporize at ~1000-1500 K (depending on composition) once 
they get close enough for temperatures to rise this high. But even before reaching this 
sublimation temperature, the collisional lifetime may drop below the orbital timescale and 
dust will quickly be smashed into pieces smaller than the blowout limit. 

This timescale (£ co ii) for inter-grain collisions to produce fragments smaller than the 
blowout limit is very short at these orbital radii and number densities. For a completely 
coupled narrow ring of width 2ae and height 2ai where all grains have the same size r gra i n , 
semimajor axis a, eccentricity e, and inclination i, t co \\ can be approximated by the inverse 
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of the collision frequency, nav, where n is the number density of particles, a is the cross- 
sectional area of one particle, and v is the relative speed between particles. This gives the 
scaling: 



^coll 



0.3 yr — 
\0.1 



P 



3 g cm 



-3 



10- b M ff 



M, 



grain 



r 



gram 



10 fimj V0.25 AU 



3.5 



(2) 



where p is the density of a dust grain, M grain is the total mass of the dust grains, and the 
approximation has been made that e ~ i (where i is in radians). With t co u = 0.3 yr for 
our typical warm systems, a steady state requires replenishing ~10~ 6 M e (or more) of dust 
every year, consuming >1000 M e of material in 1 Gyr, which is unreasonably large. 

Unrealistic mass-consumption rates implied by such steady-state models are worrying, 
and have been noted for several main-sequence, warm debris disk systems. The fractional 
luminosities L dust /L* we measure for Kepler systems with excesses also greatly exceed the 



maximum predicted by the steady-state collisional model of IWyatt et al.l (120071 ) given these 
stars' mature ages. The common logical conclusion is that dust masses at this level have 
to be transient. Several theories have been proposed to explai n such systems, inclu ding 
comet swarms, a "sup er comet " , or a recent collision (HD 69830; iBeichman et al.l 120051 ) . or 



a late-heavy bombardment-style dynamical instability forcing many icy bodies onto highly 



Wvatt et al. 


2007; 


Lisse et al. 


2012) 



Another possible explanation for these massive debris disks may be that these stars are 
all young (<100 Myr) and still in the process of clearing their p rotoplanetary disks. This 
hypothesis was offered for three of the seven stars presented by IWyatt et al.l (120071 ) with 
disks far too massive to be maintained by a steady-state collisional cascade. We inspected 
Palomar Sky Survey images and found no evidence that any of our sample stars are associated 
with the open clusters in the Kepler field, reducing the probability that this explains these 
systems. 

A further option is that the dust may be impact ejecta from asteroids or comets recently 
striking the extant exopla nets, similar to the production of giant planet rings in our Solar 
System (IBurns et al.lll999l ); however this is challenging given the super-Earths' large escape 
speeds, even if impact speeds are also high. 

What is clear is that these dust grains cannot be co-located with a planetesimal belt 
that is producing them. A massive asteroid belt at 0.2 A U or closer is an unsatisfactory dust 
source because t co n for the planetesimals is 1 Gyr (IWyatt et al.l 120071 ; iMoldovan et al. 



20101 ) and the source would now be gone. Dust also cannot simply spiral in to these close 
distances from a more distant source in steady-state. The high optical depths that we 
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calculate must be present in these systems (rangin g from r ~ 0. 01 — 0.5) are several orders 



of magnitude higher than the maximum allowed ( IWyattl 120051 ) in a steady-state inspiral 
scenario, meaning that the dust present in these quantities would destroy itself collisionally 
before PR drag has a chance to move it any significant distance. 

Given the large dust masses and the small fraction of stars that possess these excesses, 
episodic dust production mechanisms provide the most attractive explanation. Dust could 
be produced by a catastrophic collision between asteroids located further from the star, and 
we are seeing the dust as it spirals in toward the star. However, the probability of occurance 
required by this theory may not be high enough to explain the observed fraction of such 
systems. Dynamical instabilities forcing many bodies from the outer reaches of these solar 
systems onto eccentric orbits so they collide close to the star, perhaps similar in style to 
the Late-Heavy Bombardment in our Solar System, might be seen as an attractive theory. 
However, the fraction of systems observed to have 24 /im excesse s observed by Spitz er appears 
to be much too high given the short lifetime of such an event (IBooth et al.l 120091 ) . 



Discussion 



We find IR excesses around eight stars that also host planet candidates. The masses and 
luminosities of these excesses are all too large to be explained by a steady-state collisional 
cascade. Given that the WISE field coverage was not uniform in depth, and so many of these 
relatively faint stars suffered from contamination by nearby bright stars, it is challenging to 
say anything robust about the fraction of these planet-hosting stars that have excesses in our 
sample. Out of the 186 stars with W3 data of SNR>3, five have significant >5 5 excesses, or 
~3%. At first glance, this f r action appears sim ilar to the 2% of solar-type stars with warm 
excesses flWvatt et al.l 120071 : lLawler et al.l 120091 ) . 



However, the fractional luminosities of these Kepler systems w ith excesses are much 
higher than most known debris disks (Figure HI and iMoor et al.ll2006l ). Because of the larger 
distances and thus lower apparent brightness of these Kepler systems as compared to most 
known debris disk systems, only the brightest debris disks would be detected by WISE. 
Presumably the eight systems presented here are the very highest luminosity debris disks 
in the sample, and many more fainter disks are present around these planet-hosting stars. 
This would indicate that the fraction of stars possessing debris disks and Kepler-detected 
exoplanets is actually much higher than the fraction of field stars with debris disks. 



Surveys with Spitzer (IBryden et al.ll2009l ; iDodson- Robinson et al.ll201ll ) have found that 
hosting a radial velocity-detected exoplanet does not make a star more likely to host a debris 
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disk, and a recent WISE survey of exoplanets discovered by ground-based transit searches 
( IKrivov et al.ll201ll ) has found that the fraction of these s tars (~2%) to host warm debris 



disks (excesses in W3 and W4) is similar to that found by lWyatt et all (120071 ). 



The Kepler systems studied here are different from other disk-bearing exoplanet systems, 
whose planets are primarily discovered by the radial- velocity (RV) technique, and by ground- 
based transit searches. Exoplanets from the latter two techniques are mostly large and very 
close to their stars, with an average mas s of ~2 Jupiter masses and <20% known to be in 
multiplanet systems (IWright et al.ll201fl ). The recently announced 1200 Kepler exoplanet 
candidates ( IBorucki et al. are in a different regime, with about 1/3 in multiplanet 

systems, and having many super-Earths (typically 2-6 Earth radii) well inside 1 AU. Among 
our excess systems, all but two of the planets are Neptune mass or smaller; only one is 
massive enough to be a hot Jupiter. The resulting correlation we find between a star hosting 
both small planets and a debris disk inside a few AU agrees with th e planetary formation 
and dynamical evolution models presented by iRaymond et al.l (120121 ) . This makes sense in 
a picture where migrating Jupiters clean out the terrestrial systems they pass through. 

Explaining the large dust masses required to produce these excesses remains a huge 
challenge. Steady-state theories appear to fail completely, as the mass consumption rates 
required would involve destroying thousands of Earth masses of planetary material over the 
lifetime of these stars. Episodic events such as catastrophic collisions between bodies may 
provide an explanation, but these events must last long enough and occur frequently enough 
to result in ~3% of solar- type stars in this state at any given moment. Much theoretical 
work needs to be done here. 

The upcoming WISE full-sky data release will access additional systems, decrease the 
error bars on existing measurements, and provide longer-wavelength measurements for some 
of these systems to see if there is dust at larger orbital distances. Spectral data on these 
candidates would be invaluable, but the host stars are faint. Future Kepler data releases will 
contain planets on more distant orbits, perhaps yielding planets in closer proximity to these 
new debris disks. 
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Table 3. WISE Data for Debris Disk Candidates 



KOI 


Wl flux 


W2 flux 


W3 flux 


W4 flux 




[mJy] 


[mJy] 


[mJy] 


[mJy] 


469 


1.49 ± 0.05 


0.83 ± 0.03 


0.49 ± 0.07 


< 1.2 


559 


2.31 ± 0.07 


1.28 ± 0.04 


0.60 ± 0.08 


< 2.5 


871 


1.14 ± 0.04 


0.63 ± 0.02 


0.39 ± 0.08 


< 2.5 


904 


1.94 ± 0.06 


1.03 ± 0.04 


0.22 ± 0.07 a 


< 1.6 


943 


0.97 ± 0.03 


0.51 ± 0.02 


0.36 ± 0.07 


1.5 ± 0.5 a 


1020 


8.85 ± 0.27 


4.78 ± 0.15 


0.95 ± 0.10 


2.0 ± 0.5 


1099 


1.05 ± 0.04 


0.69 ± 0.03 


0.69 ± 0.09 


1.4 ± 0.5 a 


1564 


1.11 ± 0.04 


0.59 ± 0.02 


0.29 ± 0.07 


2.0 ± 0.6 



Note. — Fluxes and flux uncertainties are approximate, estimated from the given WISE 
magnitudes and magnitude errors. 



a SNR<3 



